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ABSTRACT 


The  integral  solution  of  the  sound  field  is  derived  from 
a  laminated/  interbedded  liquid-solid  half  space.  The  last 
layer  of  infinite  thickness  is  solid.  The  integral  over  wave 
number  is  conveniently  transformed  into  the  complex  wave  num¬ 
ber  plane  yielding  a  sum  of  normal  modes  of  propagation  plus 
the  sum  of  integrals  along  branch  cuts.  Waves  corresponding 
to  the  normal  modes  predominate  at  ranges  beyond  several  times 
the  water  depth  and  are  considered  in  detail.  Sample  numerical 
computations  are  presented  for  propagation  of  Rayleigh  waves 
in  shallow  water  on  the  Arctic  continental  shelf. 
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Figure  1 

Figure  2 

Figure  3 
Figure  4 

Figure  5 

Figure  6 

Figure  7 


Figure  8 


FIGURES 


Layered  model  showing  location  of  point  source  in 
a  liquid  layer. 

Layered  model  showing  location  of  detector  in  a 
liquid  layer. 

Contours  of  integration. 

Dispersion  curves  and  the  excitation  functions  of 
pressure  for  the  fundamental  Rayleigh  mode.  Com¬ 
putations  for  Model  A  of  Table  1. 

Ratio  H/e  of  horizontal  to  vertical  particle  motion. 
Computations  for  Model  A.  Positive  ratio  corres¬ 
ponds  to  retrograde  motion  and  negative  ratio  to 
prograde  motion. 

.  Dispersion  of  flexural  waves.  Computations  for 
Model  A. 

Ratio  H/m  of  horizontal  to  vertical  particle  motion 
in  the  ice  sheet  of  the  fundamental  Rayleigh  mode 
and  the  flexural  mode.  Computations  for  Model  A. 
Positive  ratio  corresponds  to  retrograde  motion  and 
negative  ratio  to  prograde  motion. 

Range  dependence  of  the  Rayleigh  moq3e  for  parameters 
of  Model  A  in  Tables  1  and  2. 


Figure 


9. 


Typical  ray  paths  for  a  source  100  m  deep  in  the  cen¬ 
tral  Arctic  Ocean.  Forty  rays  computed  at  1-degree 
intervals  at  the  source.  Angles  at  the  source  go 
from  20  degrees  above  to  20  degrees  below  the  horizon¬ 
tal. 


Figure  10.  Dispersion  and  depth  independent  excitation  for  first 

two  normal  modes  computed  for  Model  2  but  in  4  km  of 
water. 


Figure  11.  Two  models  for  computations.  Parameters  of  ice  layer 

given  in  Table  1.  Computations  for  Model  1  by  W.K.B. 
method;  those  for  Model  2  by  exact  theory. 


Figure  12.  Variation  of  pressure  and  vertical  particle  velocity 

with  depth  for  the  first  mode  at  several  frequencies. 
Computations  for  Model  2  in  4  km  of  water. 

Figure  13.  Ratio  of  horizontal  to  vertical  particle  motion  (H/z) 

at  the  ice  surface  as  a  function  of  frequency  for  hy¬ 
droacoustic  waves  of  the  first  normal  mode  and  flex¬ 
ural  waves.  Models  1  and  2  are  shown  in  Figure  10. 
Positive  ratio  corresponds  to  retrograde  motion  and 
negative  ratio  to  prograde  motion. 

Figure  14.  Variation  of  ratio  H/-&  with  depth  in  the  ice.  Positive 

ratio  corresponds  to  retrograde  motion  and  negative  ra¬ 
tio  to  prograde  motion.  Models  1  and  2  are  shown  in 
Figure  10. 
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compressional-wave  velocity  in  the  m-th  layer 

shear-wave  velocity  in  the  m-th  layer 

thickness  of  the  m-th  layer 

vertical  coordinate 

range  between  source  and  detector 

time 

angular  frequency 
phase  velocity 
wave  number 

imaginary  part  of  a  complex  number 

velocity  potential  in  the  m-th  layer 

normal  stress  in  the  m-th  layer  parallel  to  e  axis 

tangential  stress  in  the  m-th  layer 

pressure  in  the  m-th  liquid  layer 

horizontal  particle  displacement  in  the  m-th  layer  in  the 
radial  direction 


vertical  particle  displacement  in  the  m-th  layer 
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) 


horizontal  particle  velocity  in  the  m-th  layer  in  the 
radial  direction 


)OLTrw\  ^  vertical  particle  velocity  in  the  m-th  layer 
density  in  the  m-th  layer 
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Bessel  function  of 
Bessel  function 
Hankel  function  of 


«  Hankel  function  of 
Bessel  function  of 
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order  0 
of  order  0 

the  first  kind  of  order  0 
the  second  kind  of  order  0 
order  1 
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INTRODUCTION 


In  a  previous  investigation  (Kutschale,  1970)  we  derived  by  ma¬ 
trix  methods  of  Thomson  (1950),  Haskell  (1953),  Dorman  (1962),  and 
Harkrider  (1964)  the  integral  solution  of  the  wave  equation  for 
point  sources  of  harmonic  waves  in  a  liquid  layer  of  a  multilayered 
half  space  of  interbedded  liquid  and  solid  layers.  The  solution  is 
in  a  form  convenient  for  computation  on  a  high-speed  digital  compu¬ 
ter.  However,  only  the  cose  was  considered  of  a  liquid  half  space 
underlying  the  stack  of  layers.  Many  experiments  on  shallow-water 
propagation  of  low-frequency  waves  (Pekeris,  1948,  Tolstoy,  1958, 
Hunkins  and  Kutschale,  1961)  have  shown  that  this  representation  of 
the  sediments  by  a  layered  liquid  bounded  below  by  a  liquid  half 
space  is  often  useful  and  valid,  but  in  some  experiments  on  shallow- 
water  propagation  in  the  Arctic  Ocean  it  has  been  necessary  to  con¬ 
sider  the  rigidity  of  the  sediments  to  take  account  of  propagation 
of  Rayleigh  waves  in  the  layered  system.  The  ice  sheet  is  represent¬ 
ed  by  a  solid  layer  at  the  surface. 

In  the  present  report  we  extend  the  integral  solutions  of  pres¬ 
sure  and  particle  displacements  derived  from  a  layered  medium  bound¬ 
ed  below  by  a  liquid  half-space  to  a  layered  medium  bounded  below 
by  a  solid  half-space.  The  integral  solutions  of  pressure  or  particl 
displacements  are  evaluated  as  a  sum  of  normal  modes  of  propagation 
plus  the  sum  of  branch  lire  integrals  in  the  complex  wave  number 
plane.  Waves  corresponding  to  the  normal  modes  generally  predominate 
at  ranges  beyond  several  times  the  water  depth  and  will  be  considered 
in  detail.  The  branch  line  integrals  are  difficult  to  evaluate  numer 
ically  for  a  solid  half  space  underlying  the  stack  of  layers  and  dis¬ 
cussion  of  them  is  omitted.  The  formulas  for  pressure  and  particle 
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Fig.  1 
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displacements  of 
computer  programs 
liquid  half  space 
1970)  . 


the  normal  modes  are  incorporate 
for  the  IBM  360/91  used  for  the 
underlying  the  stack  of  layers 
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FORMAL  SOLUTION 
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The  Am  and  Ara  1  are  constants. 

Placing  the  origin  at  the  (m-l)-st  interface  we  have 


(1) 


(2) 


i 

Substituting  expressions  for  (Am  -  Am  )  and  (Am  +  Am ' )  from 


Equations  (1)  in  Equations  (2)  yields 
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since  the  vertical  particle  velocity  and  pressure  are  continuous 

across  each  interface. 

Consider  now  a  solid  layer  between  two  liquid  layers  or  at  the 
surface  of  the  laminated  half-space.  For  this  layer  we  may  write 
(Thomson,  1950;  Dorman,  1962) 


where  the  fourth  element  of  the  column  vector  is  the  tangential  stress 


which  is  zero  at  a  solid-liquid  boundary.  The  matrix  elements  are 
given  in  Appendix  A  and  may  be  derived  following  Haskell  (1953). 


Equation  (3)  yields  the  three  equations 
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from  Equation  (4)  in  Equations  (5)  and 
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or  in  matrix  notation 


which  is  the  same  as  the  matrix  relation  for  a  liquid  layer.  In 


general,  then,  for  interbedded  solids  and  liquids 
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H  to  4-  k 


are  solid  layers  in  sequence  we 


have 

where  each  layer  matrix  is  for  a  solid  layer.  The  product  matrix  for 


the  solid  layers  between  liquid  layers  may  then  be  written 
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Point  Source  of  Harmonic  Waves  in  a  Liquid  Layer 


LAo.+ Oa\ 


We  divide  the  liquid  source  layer  into  two  layers  as  shown  in 


i 

Figure  1.  At  z  =  D  the  pressure  is  continuous.  The  vertical  particle 


velocity  is  continuous  everywhere  in  the  plane  defined  by  z  =  D1  ex¬ 


cept  at  the  point  source  where  the  liquid  above  and  below  the  source 


moves  in  opposite  directions.  This  may  be  expressed  by  writing 
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For  the  liquid  source  layer 
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and  for  the  layers  above  z  =  D  1 


In  Equations  9  and  10 
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For  the  liquid  layers  in  terms  of  4  x  4  matrices 
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and  for  solid  layers  between  liquid  layers  in  terms  of  4  by  4  matrices 
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(Haskell ,  1953) 
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Applying  the  condition  that  no  sources  radiate  from  infinity 


From  Equation  9,  Equation  11  may  be  written 


(ID 


0 


: 


20 


^  \ 


21 


Let  v 


We  may  prove  following  Harkrider  (1964)  that  the  inverse  of  the 
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or  substituting  Equation  19  in  Equation  20 


and  solving  for  X 
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and  therefore  from  Equations  (15)  and  (16) 
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from  the  Equation  (14) 
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The  integral  solution  for  the  horizontal  particle  velocity 


at  the  surface  of  an  ice  sheet  is  by  Equation  4 
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is  the  ratio  of  horizontal  to  vertical  particle  velocity  or  particle 


displacement  at  the  surface. 


26 


n.-\ 


V 


If  a  hydrophone  or  vertical  particle  velocity  detector  is 


whe  re 
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Hence,  the  integral  solution  for  pressure  at  the  bottom  of  layer  Dx  is 
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In  particular  at  the  ocean  bottom 


♦  » 


The  corresponding  horizontal  particle  velocity  at  the  ocean 


bottom  is  obtained  from  the  exoression 
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Carrying  out  the  matrix  multiplication 
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explicit  expressions  for 


L,  K,  G,  N,  M  and  H  are 
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EVALUATION  OF  THE  INTEGRAL  SOLUTIONS 

The  integral  solutions  have  singularities  corresponding  to 

nk-mlA  -o 

k-  k<^ 

These  are  simple  poles  for  real  wave  number  and  the  integral  solu¬ 


tions  are  conveniently  evaluated  by  contour  integration  in  the 


complex 


■  -v 


plane.  The  contours  are  shown  in  Figure  3. 


First  the  transformation 


is  made.  In  addition  to  simple  poles  at 


ML  ~  o 

! 


branch  points  occur  at  kk  for  the  liquid  layers  and 


for  each  solid  layer. 


p  >  @p 

To  make  the  integrand  single  valued  we  must  introduce  branch 


cuts  originating  at  the  branch  po 


ints 


layer,  and 


Vt  .krc 


I  for  each  liquid 


for  each  solid  layer.  We  may  then 


apply  Cauchy's  Integral  Theorem.  We  note  that  since  the  integrands  are 


3? 


even  functions  of 


AO 


branch  line 


/Vv\  ^  ^  ^ 

integrals  corresponding  to  branch  cuts  made  to  these  branch  points 

cancel,  and  only  the  branch  cuts  for  ^  P  K  must  be  considered 

^  pn^ 


va.' 


The  Riemann  surface  has  four  sheets.  To  satisfy  the  vanishing 


of  VxT  as 


t 


approaches  infinity  we  must  remain  on  the 


sheet  where  the  real  part  of  \  .  and  (  f  \n 


is  greater  than 


zero.  For  the  contours  shown  in  Figure  3  complex  poles  have  been 


displaced  to  an  unused  Riemann  surface  (Ewing  et  al  1957,  pp  136-137). 

♦  cut 

The  integral  for  ur  ,  dropping  the  (p  term,  is 

now  transformed  in  the  upper  half-plane  to 
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Adding  Equations  21  and  22  and  noting  that  the  integrals  along  the 
infinite  arcs  such  as  l  l  KS\  k\\  are  zero  we  have  the  desired 
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result  ©r\ 

/-OO  / 
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The  branch  line  integrals  are  difficult  to  evaluate  numerically 
and  discussion  of  them  is  omitted,  but  we  note  that  the  proper 


signs  for  the  radicals  \  .  \  ^ 

°L  L(V 


must  be  made  on  each 


side  of  a  cut. 


kr, 


We  have  as  our  final  result  for  the  normal  modes  for 


and 


where  the  period  equation  (phase  velocity  dispersion) 


(KM-LM^-O 


has  been  used  to  simplify  the  expressions.  Group  velocity  disper- 

\J  . 


sion  is  given  by 


Corresponding  expressions  may  be  written  for  the  vertical  particle 
velocity  at  depth  or  the  horizontal  particle  velocity  at  the  ice 
surface  or  on  the  ocean  bottom. 
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LT  0^  ~  G-  ~ir  aa/\  (^m') 


For  programming  it  is  convenient  to  write  the  expressions  for 


vertical  particle  velocity  at  the  surface  and  pressure  in  the  form 


<L 


•37 


If  we  set 

=  T0ivcv<:\o^ 

in  Equations  23  and  24/  we  have 

'w. 


and  for  the  absolute  value  of  the  vertical  particle  velocity,  a 
quantity  conveniently  measured  in  transmission  experiments 


i  W0  \  = 


where 


38 


'  dc 


C-Cn 


Likewise  for  the  vertical  particle  displacement,  multiplying 


and  for  pressure 


or  for  solid  layers  between  liquid  layers 


41 


NUMERICAL  COMPUTATIONS 

2 

The  present  formulas  for  the  pressure  in  dynes/cm  and  the 
vertical  and  horizontal  particle  displacements  in  millimicrons 
are  incorporated  in  the  same  computer  programs  employed  pre¬ 
viously  (Kutschale,  1970).  Computations  are  made  in  two  stages. 
The  first  program,  an  extension  of  Dorman's  (1962)  PV7  disper¬ 
sion  program,  computes  phase-and  group-velocity  dispersion, 
the  variation  of  pressure  with  depth  in  the  water,  the  ratio  of 
horizontal  to  vertical  particle  motion  in  the  ice  and  on  the 
ocean  bottom,  the  excitation  function  dependent  only  on  the 
layering  of  the  medium  and  the  excitation  function  for  the  par¬ 
ticular  source  and  detector  depths.  For  the  1-th  normal  mode 
the  excitation  function  of  pressure  dependent  only  on  layering 
is  defined  by 

<2  jiiu)  A 

V  zL 


42 


and  the  excitation  function  for  the  particular  source  and  detector 


These  definitions  were  chosen  to  be  useful  at  long  ranges  where 


Corresponding  expressions  are  defined  for  particle  velocity  and 
displacement.  The  second  program  computes  and  plots  pressure  or 
particle  displacement  as  a  function  of  range. 

Sample  numerical  computations  are  presented  for  the  layered 
Model  A  of  Table  1.  Layer  thicknesses  and  compressional  wave 
velocities  are  similar  to  those  of  Hunkins  and  Kutschale  (1963) 
used  for  describing  explosive  sound  transmission  on  the  Chuckchi 
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TABLE  1 


Layer  Parameters  for  Model  A 


Layer 

Thickness, 

m 

Compre  ssional- 
Wave  Velocity, 
m/ sec 

Shear-Wave 

Velocity, 

m/sec 

Density , 
gm/cHW 

Ice 

1 

1 

3 

3500.0 

1800.0 

0.900 

Water 

67 

1440.0 

0.0 

1.025 

Sediment 

15 

1750.0 

250.0 

1.800 

Sed ime  nt 

Oo 

2300.0 

1 

i 

i 

i 

! 

1000 .0 

2.050 

TABLE  2 


Source  and  Detector  Parameters 


for  Model  A 


Source  level  at 

1  m,  dynes/cm^ 


Source  depth, 
m 


Detector  depth,  m 


EXCITATION  VELOCITY  (KM/SEC) 


45 


Fig.  4 


Shelf  north  of  Barrow,  Alaska.  The  ice  sheet  is  represented  by  a 
uniform  layer  3  m  thick.  The  assignment  of  shear-wave  velocities 
to  the  sediments  was  made  on  the  basis  of  observations  of  dispersion 
of  Rayleigh  waves  from  explosions  detonated  on  the  Chuckchi  Shelf 
(Hunkins,  personal  communication). 

Figure  4  shows  phase  and  group  velocity  dispersion  for  the 
fundamental  Rayleigh  mode  as  well  as  the  depth  independent  and 
depth  dependent  excitation  functions.  Source  and  detector  depths 
for  Model  A  are  given  in  Table  2.  The  computations  show  that  the 
ice  sheet  has  very  little  effect  on  amplitudes  in  the  water  and 
on  dispersion.  The  principal  effect  of  the  ice  sheet  is  to  change 
the  polarization  of  the  waves  near  the  surface.  If  the  ice  sheet 
is  neglected  only  vertical  particle  motion  exists  at  the  surface, 
but  when  the  ice  sheet  is  included  in  the  layered  system  the  waves 
are  elliptically  polarized  at  the  surface  and  change  their  orbital 
motion  through  the  ice  sheet.  This  effect  of  the  ice  on  particle 
motion  near  the  surface  is  also  observed  for  hydroacoustic  waves 


traveling  in  the  Arctic  SOFAR  channel  (Kutschale,  1972).  Figure  5 


RATIO  H/Z 
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FREQUENCY  (HZ) 


Fig.  5 
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shows  the  ratio  of  horizontal  to  vertical  particle  motion  of  the 
fundamental  Rayleigh  mode  on  the  ice  surface  as  well  as  on  the 
ocean  bottom.  For  comparison  Figure  5  also  shows  the  same  ratio 
of  flexural  waves  at  the  ice  surface.  Figure  6  shows  the  disper¬ 
sion  of  the  flexural  waves  for  Model  A.  Flexural  waves  are  surface 
waves  traveling  in  the  ice  sheet  with  pressure  decreasing  exponen¬ 
tially  with  depth  below  the  ice.  These  waves  are  weakly  excited 
from  sources  at  depth,  but  they  form  the  predominant  ice  vibrations 
from  pressure  ridging  and  fracturing  at  the  edges  of  the  floe  on 
which  listening  is  done.  Figure  7  shows  that  both  Rayleigh  waves 
and  flexural  waves  exhibit  prograde  elliptical  motion  on  the  sur¬ 
face  of  the  ice  and  retrograde  motion  on  the  bottom  of  the  ice. 
Since  the  vertical  motion  is  nearly  constant  in  the  ice  from  top  to 
bottom  for  both  flexural  waves  and  Rayleigh  waves,  Figure  7  essen¬ 
tially  shows  the  variation  of  horizontal  motion  with  depth.  Compu¬ 
tations  of  the  ratio  of  horizontal  to  vertical  particle  motion  in 
the  ice  at  frequencies  up  to  25  He  exhibit  a  similar  polarization 
to  that  of  Figure  7  for  both  types  of  waves.  On  the  basis  of  par¬ 


ticle  motion,  the  low  frequency  ice  vibrations  of  the  Rayleigh  mode 
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RANGE  (KM) 


may  be  classified  as  antisymmetric  vibrations  of  the  plate;  that 


is,  a  bending  of  the  plate  in  a  manner  similar  to 

waves. 

In  Figure  8  the  range  dependence  of  pressure 
for  the  fundamental  Rayleigh  mode.  Source  level. 


that  of  flexural 

at  5  He  is  plotted 
source  depth  and 


hydrophone  depth  are  given  in  Table  2  for  Model  A 
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Fig.  10 
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SUMMARY  OF  WORK  DONE  ON  BEHALF  OF  NOL 

The  present  report  completes  a  formulation  for  describing  low- 
frequency  propagation  in  the  ice  covered  Arctic  Ocean.  In  a  report 
in  preparation  numerous  numerical  computations  will  be  presented  for 
deep  and  shallow  water  propagation  to  illustrate  the  effects  of  an 
ice  layer  on  group  and  phase  velocity  dispersion,  variation  of  pres¬ 
sure  with  depth  in  the  ocean,  particle  motions  in  the  ice,  and  the 
variation  of  pressure  and  particle  motions  as  a  function  of  range 
from  the  source.  We  summarize  the  principal  features  here  in  deep 
water  for  waves  propagating  in  the  SOFAR  channel.  Sample  ray  paths 
are  shown  in  Figure  9.  Figure  10  shows  phase  and  group  velocity  dis¬ 
persion  and  the  depth  independent  excitation  function  for  the  first 
two  modes  computed  for  Model  2  of  Figure  11  but  in  water  4KM  deep. 

Tne  variation  of  pressure  and  vertical  particle  velocity  in  the  water 
for  the  first  mode  is  shown  in  Figure  12.  Parameters  for  the  ice 

sheet  are  the  same  as  in  Table  1. 

For  an  ice  layer  3  m  thick,  typical  of  the  central  Arctic  Ocean, 
the  effect  on  pressure  amplitudes  at  depth  and  dispersion  is  small, 
but  even  this  thin  ice  sheet  causes  a  large  change  in  particle  mo¬ 
tions  near  the  surface.  The  effect  is  similar  to  that  of  Rayleigh 
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RATIO  H/Z,  40  HZ 
5  0  .5 


waves  illustrated  in  the  previous  section.  Waves  are  elliptically 
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polarized  in  the  plane  of  propagation.  The  ratio  of  horizontal  to 
vertical  particle  motion  is  quite  sensitive  to  ice  thickness.  To 
speed  the  numerical  work,  a  rapid  method  was  developed  for  compu¬ 
ting  particle  in  the  ice  by  the  W.K.B.  method  (Kutschale,  1972). 
Figure  13  shows  an  example  of  the  ratio  of  horizontal  to  vertical 
particle  motion  at  the  ice  surface  for  hydroacoustic  waves  of  the 
first  normal  mode  computed  by  the  W.K.B.  method  (Model  1  of  Fig.  11) 
and  exact  theory  (Model  2  of  Fig.  11).  For  comparison  this  same 
ratio  is  shown  for  flexural  waves.  Ice  thickness  is  3  m.  A  posi¬ 
tive  ratio  corresponds  to  retrograde  elliptical  motion  and  a  nega¬ 
tive  one  to  prograde  elliptical  motion.  Flexural  waves  generated 
by  large-scale  ice  movements  at  the  boundaries  of  floes  are  the 
principal  background  noise  in  the  ice  and  it  is  important  to  com¬ 
pare  the  particle  motions  of  flexural  waves  in  the  ice  with  those 
of  the  hydroacoustic  waves.  Fig.  14  shows  the  principal  features 
at  low  frequencies  for  both  types  of  waves.  Particle  motion  is  re¬ 
trograde  elliptical  at  the  surface  and  prograde  elliptical  at  the 
bottom  of  the  ice.  Note  that  Rayleigh  waves  of  Figures  5  and  7 


also  show  this  type  of  particle  motion.  since  the  vertical  motion 


is  nearly  constant  through  the  ice  for  both  flexural  waves  and  hy¬ 
droacoustic  waves,  these  plots  essentially  show  the  variation  of 


horizontal  particle  motion  with  depth  in  the  ice.  It  is  signifi¬ 
cant  that  a  node  of  horizontal  particle  motion  often  does  not  occur 
at  the  same  depth  in  the  ice  for  both  hydroacoustic  waves  and  flex¬ 
ural  waves.  This  effect  may  be  useful  to  improve  the  signal  to 
noise  ratio  when  horizontal  component  geophones  are  used  as  listen¬ 
ing  devices. 

Future  work  will  emphasize  analysis  of  field  data  and  compari¬ 
son  of  these  data  with  theory.  Over  one  hundred  recordings  from  ex¬ 
plosions  have  been  made  in  deep  water  and  shallow  water  areas  under 
various  ice  conditions  employing  as  listening  devices  vertical  and 
horizontal  component  geophones  on  the  ice  as  well  as  hydrophones  at 
depth.  An  analysis  of  these  recordings  will  establish  the  effects 
of  variations  of  ice  thickness  and  roughness  on  acoustic  wave  veloci 
ties,  wave  amplitudes,  and  particle  motions  on  the  ice.  Guided  by 
these  results  and  the  numerical  work  additional,  experiments  will  be 
made  employing  explosive  and  harmonic  sources  to  verify  the  theoreti 
cal  predictions  on  the  variation  of  particle  motion  of  hydroacoustic 
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and  flexural  waves  as  a  function  of  depth  in  the  ice  and  to  determine 
whether  in  fact  the  signal  to  noise  ratio  can  be  enhanced  by  a  hori¬ 
zontal  component  geophone  buried  near  mid-depth  in  the  ice. 

Since  large  scale  ice  movements  are  the  principal  source  of 
background  noise  we  will  investigate  the  spectral  character  of  the 
ice  vibrations  generated  by  these  movements.  Data  have  been  record¬ 
ed  from  geophones  and  hydrophones  installed  in  an  active  pressure 
ridge  and  additional  experiments  are  planned. 
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APPENDIX  A 
MATRIX  ELEMENTS 


Matrix  elements  (see  Haskell  1953,  and  Dorman,-  1962).  Solid  layers 

CQ-Ou  -  C 

(a  ^  -  (a  =  <:  ^  rn 

CA-O,^  ~  ~  ~ ^  ( c 

(a^,--^--  -C  i>vrv^^  +- 

r^K-^oaQ^ 

t<W)w  -  (a»Ow  =  -  c  w-  cosU,  + 
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Liquid  layers  4x4  matrices 
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